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Many independent high-resolution simulations of structure formation in cold dark matter models 
show that galactic halos should have singular core profiles. This is in stark contrast with observations 
of both low- and high-surface brightness galaxies, which indicate that the dark matter has almost 
constant density in the central parts of halos. Basically there are three possible avenues to a solution 
to the problem, which we discuss in turn. Observations of halo profiles could be more uncertain 
than previously thought, and higher resolution observations could reveal that spirals do have a 
singular core feature. The highest resolution simulations do not include a baryonic component, 
and it is conceivable that violent star formation processes and similar phenomena can destroy the 
singular dark matter core and lead to an almost constant density core profile. Finally, we discuss 
in more detail the intriguing possibility that the discrepancy hints at some new exotic physics of 
the dark matter. Warm dark matter and self-interacting dark matter are two of the most promising 
candidates. 



I. INTRODUCTION 

Dark matter seems to be a necessary ingredient for structure formation in the universe |lj . The standard hierarchical 
clustering model where the dark matter is in the form of very non-relativistic, collisionless particles has been very 
successful in explaining structure from galactic scales to the largest scale observable, that of the cosmic microwave 
background radiation (CMBR) In the 1980s it was realised that the flat, pure CDM model produces too much 

small scale structure if it is normalized on large scales js). Some modification of the model is needed, and several 
possible avenues for this exist. The most popular model initially was mixed dark matter where some of the dark 
matter is made up of light particles with large free streaming length. The prime candidate for these light particles is 
a neutrino with a mass of a few eV. However, recent results from Super-Kamiokande Q and other neutrino oscillation 
experiments suggest that the mass of the most massive neutrino is of the order 0.1 eV, far too small a mass to 
be cosmologically significant. Thus, the mixed dark matter model is currently disfavoured although perhaps not 
completely excluded. However, data from type la supernovae at high redshift suggests that the energy density in the 
universe is dominated by a cosmological constant ||J7|. This also dampens small scale structure formation and can 
bring the CDM model into agreement with observations without the need for hot dark matter. 

Recently, however, a new problem for the CDM model has surfaced, having to do with structure on galactic scales 
and below. Observations of spiral galaxies, both high and low surface brightness, suggest that the dark matter in 
spirals is in the form of a halo with almost constant central density ||s|-|l2|. To a reasonable approximation it can be 
described by a universal profile of the form jl2| 

P( r ) - Poi — ; T7~2 i — 2T -V-rt Pa, (1) 

where the central density, po, and the core radius, ro, depend on the specific system. However, numerical simulations 
of halo formation consistently show much more singular core profiles. The first large scale study of CDM halo profiles, 



with sufficient resolution to resolve galaxy sized halos, was that of Navarro, Frenk and White |13|, who found that 
the simulated CDM halos also follow a universal profile 

p(r) ^ i , w. c , c ] s 2 ~>r^0 PcSc—, (2) 
(r/r s )(l + r/r s y r 

known as the NFW profile. p c 8 c and r s are again parameters that are to be fitted to individual systems. Close to 
the center this profile approaches p oc r -1 , in strong disagreement with observations. More recent simulations with 
much higher resolution find an even steeper profile, approaching p oc r~ 3 / 2 p|-|l7|]. Furthermore, it has been claimed 
that these CDM simulations are now of sufficiently high resolution that the results for the central halo profile have 
converged to the infinite resolution limit. 

In addition to this problem with the halo density profiles, observations show that a galaxy like the Milky Way 
contains almost an order of magnitude fewer small satellite galaxies than what is found in CDM simulations [|bg|]l9]l . 
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II. POSSIBLE SOLUTIONS 



This discrepancy between CDM simulations and observations of galaxies is so serious that it needs to be investigated 
more closely, because it could indicate a fundamental flaw in our understanding of structure formation and the nature 
of dark matter. There are basically three possible ways to remedy the discrepancy, which we shall discuss in turn. 

A. Observational solutions 

Until very recently, relatively little effort has been devoted to measuring reliable rotation curves in the inner parts 
of low surface brightness (LSB) galaxies, the reason being that the outer parts of the rotation curves are completely 
dominated by dark matter, and thus hold important information about the total dark matter mass in the halo. Even 
if there is a singular core in the halo, its integrated mass is only a small fraction of the total halo mass. In that case, 
high angular resolution in the central part of the galaxy is not an important issue. Therefore observations of galactic 
rotation curves have primarily been made using observations of the 21cm neutral hydrogen hyperfine-structure line 
with radio telescopes. These observations have the advantage that neutral hydrogen systems extend further out in 
most galaxies than the visible light, meaning that the halo structure can be mapped out to very large radii. If the total 
halo mass is the parameter to be determined then this technique is indeed the best possible. But because of the poor 
angular resolution of radio observations there are indeed significant effects of beam-smearing in the measurements of 
the inner rotation curves. Furthermore there is usually not much neutral hydrogen present within the central kpc of 
a galaxy. 

Recently it was pointed out by van den Bosch and Swaters |20| that taking beam smearing into account increases 
the uncertainty in rotation curve measurements by a significant amount. In fact it seems that some LSB galaxies are 
consistent with having a cusp-like feature in the central halo, whereas others are still best fitted by a constant density 
core. At present it is not clear how important this beam-smearing effect really is, and other very recent studies reach 
the opposite conclusion, that galaxy halos are best fitted by constant density cores |ll[| . 

Instead of using radio observations to probe the inner parts of rotation curves it may be more feasible to use 
observations of HII regions instead of neutral hydrogen. These observations can be made with optical telescopes and 
in principle have very high angular resolution. HII is usually also very abundant in the central parts of galaxies. One 
very recent study of HII rotation curves by Swaters, Madore and Trewhella Q indicates that there could be cusps 
present, but on the other hand they also do not rule out the constant density cores. Clearly, better observational 
data from the inner parts of galaxies are needed, as well as an increased understanding of the possible uncertainties 
involved in these observations. 

It is more difficult to imagine an observational solution to the substructure problem. If small satellite galaxies were 
present in galactic halos in the numbers predicted by CDM simulations they would doubtless have been detected. 
However, that conclusion assumes that their mass-to-light ratio is not too different from ordinary galaxies. As will 
be discussed in the next section, it has been proposed that small satellite galaxies are not seen, simply because they 
have extremely large mass-to-light ratios. 

B. Astrophysical Solutions 

CDM simulations with very high resolution (sufficient to resolve cores of galactic halos) at present do not include 
baryons jl^Jl^| . The reason is both that the computational demands for this type of simulation are very high, and 
that there is no really reliable way of including star formation, supernovae and similar small scale features in the 
baryonic component. Without such high resolution simulations it is difficult to completely exclude the possibility that 
a proper treatment of the baryons would cure the discrepancy between simulations and observations. 

Several scenarios have in fact been proposed that try to remedy the problem by invoking baryonic features. As 
was mentioned in the last section, small satellite galaxies would not have been detected if they contain very few 
stars. Bullock, Kravtsov and Weinberg J2^] have suggested that when the universe was reionised at a redshift z > 5, 
accretion of baryons onto existing low mass dark matter halos was almost completely halted. The presence of a strong 
ionising radiation field stops formation of galaxies with velocity dispersion smaller than ~ 30 km/s, almost exactly 
the scale below which the discrepancy between the predicted and observed number of systems becomes apparent. If 
this explanation is true, then the prediction is that galactic halos contain a very large number of dark matter clumps 
with very few stars in them. Such clumps could possibly have disastrous effect on disk galaxies. It is well known that 
thin disks are unstable to heating by the frequent passage through the disk of massive objects, whereas thick disks are 
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more robust. At present it is not clear how the distribution of dark matter clumps predicted by CDM would affect 
disk galaxies, and better simulations of this effect would be highly desirable. 

It has also been proposed that tidal interactions between baryons and dark matter in the central halo can erase the 
cusp-like feature seen in pure CDM simulations |23j . In this scenario a singular CDM halo forms at high redshift and 
starts accreting a baryonic disk. At some point a massive burst of star formation and subsequent supernova events 
occurs, initiating a galactic wind that can eventually carry off almost all gas present in the disk. If not too much star 
formation has occurred prior to the wind formation, the total baryonic mass carried off could be close to the total 
mass of the disk. This dramatically changes the gravitational potential felt by the dark matter and the subsequent 
violent relaxation can erase the singular core. After this event, the baryonic disk which is observed in present day 
spirals is slowly accreted at lower redshift. 

Even though this scenario seems plausible there are a number of serious questions related to it. Firstly, preliminary 
results of galactic winds driven by star formation show that almost all metals can be easily expelled from the galaxy 
because of their high opacity |24|,^5|]. This fits well with the fact that metals are observed in abundance in the 
intergalactic medium at high redshift. However, it seems impossible to expel a large fraction of the hydrogen, which 
makes up most of the mass |24[]25| 1 . Thus it seems unlikely that the required amount of baryons can be expelled from 
the galaxy by this initial burst of star formation. Secondly, the hypothesis that most of the baryons remain in diffuse 
gas until they are expelled is neither supported nor ruled out by present observations. In order to determine if this is 
indeed the better understanding of star formation is needed. 

Thus, the final conclusion is that a scenario where the problem is resolved by properly taking into account inter- 
actions between baryons and the dark matter cannot be ruled out at present. However, explaining the cusp problem 
requires some fairly strong assumptions about the behaviour of the baryonic component at high redshift. 

C. Possible particle physics solutions 

If the astrophysical or observational solutions turn out not to work, then the only option we have left is to consider 
solutions from particle physics. Roughly this type of solution can be divided up into two categories, those that derive 
from physics at the epoch of initial power spectrum formation and those which come from changing the physical 
properties of the dark matter particles. 

The first type of solution is in a sense the simplest because it does not affect physics during structure formation, 
only the initial conditions. The initial condition in this regard is the power spectrum of fluctuations from inflation. In 
order to solve the substructure problem we need to dramatically reduce fluctuations on scales smaller than a galaxy, 
or 

k > 5 f M 
1 Mpc" 1 ~ {W^Mq 

This could for instance be the case in inflationary models with broken scale invariance, as was proposed by 
Kamionkowski and Liddle |2(|. These authors showed that the substructure problem could indeed be solved by 
this method. However, the cusp problem is more serious. It has been shown in N-body CDM simulations, with an 
artificial power spectrum cut-off introduced, that dark matter halos also have singular cores 

@- The reason for this 

is simply that in any cold, collisionless collapse there is a large amount of material with very low entropy (tempera- 
ture). This material eventually ends up in the centers of gravitational wells, forming the singular cores observed in 
simulations. Thus, the conclusion is that introducing a power spectrum cut-off does not remedy the problem, it is a 
generic feature in any collisionless CDM structure formation scenario, regardless of the initial conditions. 

This leads us to consider the other possible solution which has to do with the nature of the dark matter itself. A 
very interesting possibility was proposed by Spergel and Steinhardt [^7j (see also |2^-Q ) , namely that the cold dark 
matter could possess relatively strong self-interactions. If there are such self-interactions they have significant effect 
on halo formation. For example substructure clumps will be evaporated by interactions with the smooth dark matter 
halo background, and could decrease the number of substructure halos to the observed level. 

The maximum entropy state of any self-gravitating system is the singular isothermal sphere, regardless of possible 
self-interactions. However, if the interaction cross-section has the right magnitude, collapse to a singular halo could 
be prevented for a sufficiently long time that present day halos would not show singular cores. It has been shown in 
numerical simulations p5|-p7t that a cross section higher than 

- > 1CT 23 - 1(T 22 cm 2 /GeV (4) 
m 

leads to quick core-collapse and formation of a singular core. On the other hand, cross sections smaller than 
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— < 10~ 25 cm 2 /GeV (5) 
m 

will not have an observable effect on structure formation. If the cross section is fine-tuned to a value 

- > 10- 24 - 10- 23 cm 2 /GeV, (6) 

TO 

the initial expansion of the halo core prior to collapse lasts more than a Hubble time, so that present day halos could 
fit observations (see A. Burkert in this volume for a more thorough discussion of this model). At present it is not 
entirely clear whether the scattering cross-sections needed to explain the substructure problem are the the same as 
are needed to explain the cusp problem. If this turns out not to be the case, then self-interacting cold dark matter is 
certainly disfavoured. 

Another possibility which has received a lot of attention recently is that dark matter is not completely cold, but 
has some thermal motion in the early universe pq-ji^]. If the dark matter has thermal motion, any fluctuations are 
erased on scales smaller than the free-streaming length |ij 



^streaming 



1 Mpc" 1 (TO/lkeV) 3/4 , (7) 



for particles which are in a relativistically decoupled thermal distribution. In the old hot dark matter model of 
structure formation, where the dark matter was in the form of neutrinos with a mass of ~ 30 eV, the free-streaming 
length would be of order ^streaming — 0.05 Mpc -1 , roughly the size of a large cluster. This is clearly too large to fit 
large scale observations, but if the mass of the dark matter particle is ~ 1 keV, the free-streaming length is of the 
order ^streaming — 1 — 2 Mpc -1 , just what is needed to explain the substructure problem. At first sight this "warm" 
dark matter (WDM) model |38|-|45) would seem to suffer from the same problem as those with a power spectrum 



cut-off, namely that low entropy material clusters in the central parts of halos and produce singular cores. However, 
there is one vital difference between this model and CDM with a power spectrum cut-off. The warm dark matter has 
finite temperature, which implies that the phase-space distribution occupies a finite volume in momentum space. On 
the other hand, the CDM distribution is by definition a 5-function in momentum space. A useful quantity describing 
this is the average phase space density [ p8| , ^9| 

Q = P f < oo for WDM /on 
W - ( w 2^3/2 1 = 00 for CDM ' 1 ' 

In an adiabatic collapse to a lowered isothermal sphere, the central halo density can be crudely estimated as ]3§| , |3l| 

_ na^\3/2 / < oo for WDM , Q s 
P^Q{Ze)> | soo forCDM . (9) 

Thus, WDM halos are expected to have non-singular cores. However, the question is whether the same WDM particle 
mass can solve both the substructure problem and the cusp problem. From studies of galactic substructure, as well 
as the Lyman-a forest 0] it has been concluded that towdm si 750 eV in order to fit observations. 

To date no N-body simulation with sufficient resolution to tackle the WDM effect on halo cores has been performed. 
We are therefore for the moment forced to rely on crude estimates of the effect, such as Eq. Mj. Hogan and Dalcanton 
9[ have estimated that the WDM mass should be Towdm — 200 — 300 eV in order to produce the correct core 
ii for low surface brightness galaxies. 
It therefore seems that the collisionless WDM model has difficulty in explaining both the substructure and the cusp 
problem (however, see also Ref. pq). However, it should not be immediately dismissed because whereas the estimate 
of the mass needed to reproduce substructure is probably quite accurate, the estimate of the mass needed to resolve 
the cusp problem is at best a crude guess. 

If it turns out that collisionless WDM does not work then WDM with self- interactions might jlTj . Even relatively 
weak self-interactions will keep the WDM distribution in pressure equilibrium until T <C to in the early universe. In 
that case the WDM particles cannot free-stream and perturbation on small scales are not damped. However, the 
model does not resemble CDM either because the pressure equilibrium prevent fluctuations smaller than the Jeans 
mass from growing |4g|-|5"o|| . Since the Jeans length is almost a factor of 2 smaller than the free streaming scale |47fl 

kj ~ 1.7 to 3 c ^ Mpc- 1 (10) 

^stream — 1-1 mj^y MpC~\ (11) 



racl 



we find that the mass needed to explain the substructure problem in this model is only to ~ 300 — 400 eV, which is 
quite close to the mass needed to explain finite cores. 
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III. DISCUSSION 



We have reviewed the nature of the apparent discrepancy between CDM simulations of galactic halos and observa- 
tions of spiral galaxies. In addition we have discussed which possible solutions there are to this problem. The simplest 
possibility is that either observations or simulations are in error, but presently it is not at all clear whether or not 
this is indeed the case. If it turns out that the discrepancy persists even with better simulations and observations, 
the only other solution is some additional particle physics. Merely changing the initial power spectrum from inflation 
does not work, so the only option is that the dark matter is not standard CDM. Both self-interacting CDM and warm 
dark matter arc intriguing possibilities, although at the moment these models have not been sufficiently investigated 
to determine whether they work better than CDM. In any case, the CDM problem on small scales has stirred a lot of 
interest recently and is likely to do so in the future as well. 
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